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Abstract 

The Lieb-Oxford bound is a constraint upon approximate exchange-correlation functionals. We 
explore a non-empirical tightening of that bound in both universal and electron-number-dependent 
form. The test functional is PBE. Regarding both atomization energies (slightly worsened) and 
bond lengths (slightly bettered), we find the PBE functional to be remarkably insensitive to the 
value of the Lieb-Oxford bound. This both rationalizes the use of the original Lieb-Oxford constant 
in PBE and suggests that enhancement factors more sensitive to sharpened constraints await 
discovery. 
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I. BACKGROUND 



Construction of approximate exchange-correlation (XC) functionals in DFT without re- 
hance on empirical data is an important task, both conceptually and practically. Perhaps 
the most widely used constraint-based approximate XC functional today is the extremely 
popular Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA) 3|. One 
of the constraints on which the PBE GGA X functional (and some others also) is based is 
the Lieb-Oxford bound [5|]. In the DFT literature this bound commonly is expressed as 



ELDA^n] 

where 



< \lo (1) 



E^^^M = -U^-\" j dhnn-). (2) 



The LO value for the constant is 

\lo = 2.273 (3) 

The possibility of tightening this bound has been the subject of recurrent interest in DFT. A 
slightly tighter value Xch = 2.215 was found by Chan and Handy 6]. Vela later reported 
that using a spatially varying implementation of the LO bound which always is tighter than 
Xlo improved the results for a test set of light inorganic and organic omolecules calculated 
using constraint-based GGAs. 

Shortly thereafter and independently, two of us (OC hereafter) 8|, |9| gave numerical evi- 
dence from exact and near-exact calculations on atoms, small molecules, and model systems 
that the true bound is much tighter. That analysis proceeded by defining the functional 

AM 

with both numerator and denominator evaluated on the actual density of each system. In 
general, this functional cannot be evaluated exactly, because neither E^c nor the density 
is known exactly. However, it can be evaluated to high accuracy for systems for which 
near-exact XC energies and system densities are known from configuration interaction or 
quantum monte carlo calculations. The result js, ^ is that real systems typically have 
AM £ (1-1 • • • 1-3). The higher end of the interval typically corresponds to more rarefied, 
diffuse density distributions, while the lower end corresponds to more compact densities. 



2 



Values above 1.3 were only found for extreme low-density limits of model Hamiltonians: 
the k —>■ limit of Hooke's atom has X[k — > 0] = 1.489, and the ^ oo limit of the 
homogeneous electron gas has X[rs — oo] = 1.9555 =: Xheg- 

On the basis of these results, OC conjectured j^, 9] that further tightening of the LO 
bound, beyond that obtained by Chan and Handy, can be achieved, and suggested that for 
real systems (excluding unphysical limits of model Hamiltonians) Xoci = 1-35 may provide 
the tightest upper limit, whereas for arbitrary systems Xoc2 = 2.00 ~ Xheg is the upper 
limit. 

OC also speculated that system-specific upper limits could be found, thereby providing 
upper limits for all systems sharing some common properties. Earlier there was other ev- 
idence for system-specific limits. Novikov et al. [lO[] used a reduced k parameter (defined 
below) in the PBE XC functional to some benefit. This reduction (see our discussion below) 
is equivalent to a reduced LO bound. The numerical rationalization for this was published 
somewhat later by Peltzer y Blanca et al. [llj. Translating to effective values of A, broadly 
they found that 3d metals do better with A ~ Alo? 4d metals benefit from A ~ 1.81 —>■ 1.94, 
and 5d metals benefit from A ~ 1.69 — * 1.84. The notable exception was Fe, where the 
effective A was 2.8, an illustration of the fact that all the limitations of a specified XC form 
cannot be corrected by a single parameter fix. (Recently there also has been study of re- 
duced K in the PBE functional but the reduction is done in such a way as to respect the 
original Lieb-Oxford bound |l2|, hence is not directly related to the issue at hand.) 

Other than this one empirical example, the available data did not allow any general 
characterization of A-value classes. Here we propose and explore a generally applicable, 
entirely non-empirical way to characterize classes of systems with a common maximum 
value of X[n]. This characterization is based on a rarely mentioned part of the original 
Lieb-Oxford paper, in which they show that tighter estimates of the upper limit on X[n] 
can be achieved by restricting the A functional to densities which integrate to a specified 
particle number A^. We therefore introduce the function X{N), which for a given value of 
N provides a universal upper limit upon X[n] valid for all systems such that J d^rn{r) = N. 
The maximum value of A(A^), attained for cxd, is the value Xlo used in common density 
functionals. The function A(A^) assigns to each class of systems of common particle number 
an upper limit A(A^) < Xlo- 

In construction of constraint-based functionals, the fact that the upper limit can be 
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tightened globally (from Alo to Xch and perhaps on to A0C2) or in a system-specific way 
(e.g., using A(A^)) has not been taken into account, and the consequences of a replacement 
of Xlo by one of the lower values in currently popular functionals are unknown. We study 
some of those consequences here. 



II. CONSTRUCTION OF A PARTICLE-NUMBER DEPENDENT BOUND 



To explore the system-specific bound provided by the function A(A^) requires facing the 
problem that, while Lieb and Oxford proved the existence of this function and deduced 
some of its properties, they did not obtain a closed analytical expression for all A^. We thus 
propose a simple approximation to A(A^), compatible with all known information on the 

nnn 

universal LO bound. The following facts are known about A(A^) [5|, |8|, |9|: 

(i) Its value at = 1 is A(A^ = 1) = 1-48 := Ai. 

(ii) Its value at = 2 is not known, but is above \min{N = 2) = 1.67. 

(iii) The function A(A^) is monotonic, i.e., A(A^ + 1) > '^(^)- 

(iv) Its value at = 00 is not known, but must be less than or equal to \max{N — >■ 
00) := Aoo- Different proposals for the value of Aoo are \lo = 2.273, \ch = 2.215, and 
^oc2 = 2.00 ~ Xheg- 

(v) The largest value of X[n] found for any system studied specifically is that for the extreme 
low-density limit of the homogeneous electron gas XHEcifs 00) = 1.9555. For real 
physical systems, X[n\ typically < 1.3. These values provide empirical lower bounds on the 
function A(A^). 

Note that standard density functionals either do not make use of the Lieb-Oxford bound 



at all (and some are known to violate it [13|, I14J . I15I . Il6| ) or exploit only property (iv), 
normally with the weakest value for Aoo, namely Xlo- To construct a model for the function 
A(A^) we exploit properties (i) (value at A^ = 1), (iii) (monotonicity) and (iv) (value at 
A^ —>■ 00). We use properties (ii) (theoretical lower limit at A^ = 2) and (v) {X[n] for model 
and real systems) as consistency tests for the construction. With all this in mind, we propose 
the simple interpolation 

A(iV)=(l-^)A„ + |, (5) 

where Aoo is Xlo, ^ch or Xoc2- By construction this interpolation obeys properties (i), (iii) 
and (iv). Direct inspection shows that it also respects properties (ii) and (v). 
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FIG. 1: [Color] Interpolation function A(A^) for three different choices of Aqo, compared to known 
theoretical results and empirical data. Black crosses represent the exact value at = 1 and the 
lower bound at = 2. Values at — > cxo are indicated on the right vertical axis. The three 
continuous curves are our interpolation, Eq. (I5|), using the three alternative choices for Aqo- All 
other data represent ranges or values for selected real systems, providing empirical lower bounds. 

Figure [1] illustrates this function for the three different choices of Aoo, and compares it 
to the known value at = 1, the lower limit at A^ = 2, and some representative data for 
atoms, molecules, and the homogeneous electron gas. 

III. IMPLEMENTATION AND COMPUTATIONAL PROTOCOLS 
A. Modification of PBE GGA 



To explore these ideas we have implemented the various possible replacements of A^o in 
the PBE GGA. At the outset, we remark that, on the basis of previous experience with the 
revPBE functional jl6|, we expect that lowering of A in PBE will have a detrimental effect 



5 



on atomic total energies. (In revPBE an increase of A was shown to improve atomic total 
energies and molecular atomization energies, at the expense of worsened bond lengths.) 

Since the actual values of A[n] for physical systems are known to fall far below \loi 
and the theoretical information available from the CH numerical tightening and from the 
function A(A^) both indicate that lower values of A are appropriate, this detrimental effect 
must be considered a severe shortcoming of the GGA. An important issue of energetics, 
therefore, is whether the atomization energies are improved when tightened LO bounds are 



used in a GGA. (Further investigation is needed to see if meta-GGA functionals j4| suffer 
from the same problem, but that is beyond the scope of this study.) 

In any event, there are five possibilities for tightening, pertaining to two categories. Cat- 
egory I is a simple replacement of the constant value A lo by the alternative lower constants 
\cH or \oc2- Category II replaces the constant by the function A(A^), with the three possi- 
ble choices for Aqo- The resulting five choices are to be compared to the original choice Xlo-, 
made in the construction of PBE. 

In the PBE GGA, the LO bound is enforced pointwise through the choice of the parameter 
K in the exchange enhancement factor 



F!^^:=1 + k- ^ (6) 
1 + hs^/k ^ ' 

with the dimensionless reduced gradient given by 

.(r) = ^ (7) 

^ ^ 2(37r2)i/3 n(r)4/3 

Taking spin-polarization into account, satisfaction of Eq. ([2]) by the enhancement factor ([6]) 
is equivalent to 

F.K.]<^ = 1.804 (8) 
Since lims^oo Fx\n, s] = 1 + k, the result is 

KPBE = 0.804 (9) 

Of course, the simple choice of a different universal bound leads to 

«^(Aoo) = ^ - 1 (10) 
The equivalent modification to include the A^-dependent LO bound, Eq. ([5]), is 

K(Ar,A^) = ^i^-l (11) 



The result of considering such altered LO bounds is five variants of the PBE X functional: 

PBEMA: PBE96 X but with Aoo = Xch = 2.215. 

PBEMB: PBE96 X but with Aoo = ^oc2 = 2.00. 

PBEMC: PBE96 X but with A(A^, Aoo) and Aoo = >^lo = 2.273. 

PBEMD: PBE96 X but with A(A^, Aoo) and Aoo = ^ch = 2.215. 

PBEME: PBE96 X but with A(A^, Aoo) and Aoo = \oc2 = 2.00. 

The first two, along with the original PBE, comprise category I, the latter three, category 
II. We denote the five variants collectively as PBEMx in what follows. 



B. Protocols 



All five variants were introduced in the code deMon2k, version 2.4.2 [17], by systematic 
modification of the exchange-correlation modules. Subsequently, the implementation was 
validated by c omp arison of atomic calculations done with hard-coded modifications of the 



code soatom.f ISj. Throughout we used the full PBE correlation functional, not the deMon 
cutoff version (i.e. we used the deMon2k "PBESSF" option), both for ordinary PBE and 
PBEMx. Because deMon2k uses variational Coulomb fitting, there is a choice of density 
fitting (auxiliary) basis sets and of the method for evaluating XC matrix elements. Initially 
we used the so-called A2 density fitting basis (deMon2k option AUXIS(2)) and the option to 
do the numerical integrals for the XC quantities using the fitted (auxiliary) density (deMon2k 
"AUXIS" option). We return to these options below. 

For development of a suitable protocol (Kohn-Sham basis, fitting basis) we first studied 
the Li2 molecule in a triple-zeta-plus-polarization (TZVP) KS basis. The results are in Table 
[H Note that AE is the total atomization energy, 2 En atom — (not the cohesive energy 
per atom). Regarding the quality of the calculation, observe that for the unmodified PBE 



functional, our results are almost identical with those given by Ernzerhof and Scuseria |19| . 
AE = 20 kcal/mol, = 2.727 A. (For reference, they quote experimental values as 26 
kcal/mol and 2.673 A.) 

As would be expected from naive use of a particle-number-dependent model, the results in 
Tabled] show a clear size-inconsistency problem, signaled by the big shift in AE between the 
A^-independent models, PBE, PBEMA, PBEMB, and the A^-dependent models, PBEMC, 
PBEMD, and PBEME. The fact that there is no such shift in the Re values is a clear sign 
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that the problem is in the comparison with the isolated atom. Eq. (fTTi) illustrates the point. 
In a naive application of the A^-dependent models, the Li atom has A(3, Aoo) while the Li2 
molecule has A(6, Aoo) (with the three choices of Aqo)- The result is a separated atom limit 
of the diatomic molecule which is not the same as the isolated atom. Table [III shows the 
very substantial difference in the PBE parameter k for these two situations. 

TABLE I: Comparison of effects of various Lieb-Oxford bounds in the PBE X functional for the Li2 
molecule. See text for notation about functionals. Eatom and are total energies in Hartrees. 
/S.E is the total atomization energy in kcal/mol, is the equilibrium bond length in A. 



Functional Eatom. Eli^ AE 

PBE -7.460992748 -14.953949056 20.06 2.7236 

PBEMA -7.457436406 -14.946932892 20.01 2.7218 

PBEMB -7.441633876 -14.915640294 20.31 2.7155 

PBEMC -7.442385785 -14.937217688 32.91 2.7196 

PBEMD -7.439006265 -14.930264047 32.79 2.7181 

PBEME -7.424518016 -14.899876726 31.90 2.7131 



TABLE IL Values of the PBE X functional parameter k{N, Aoo) for the PBE and PBEMx func- 
tionals for = 3 (Li atom) and = 6 (Li2 molecule). 

Functional k(3, Aoo) '^(S, Aoo) 

PBE 0.804319 0.804319 

PBEMA 0.757967 0.757967 

PBEMB 0.587401 0.587401 

PBEMC 0.594439 0.699379 

PBEMD 0.563537 0.660752 

PBEME 0.449826 0.518614 



Table IIIII shows how a size-consistent set of parameters, here for = 6, resolves the 
problem. (For clarity, note we made the common choice throughout all these calculations 
and ignored the DFT spin-symmetry problem. Thus, the separated atoms are spin-polarized 
even though the molecule has multiplicity equal one.) Throughout this study, we used 
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TABLE III: Comparison of effects of various Lieb-Oxford bounds in the PBE X functional for the 
Li2 molecule. A^-dependent functionals done with = 6 size-consistent parameters. TZVP KS 
basis. Upper set is the A2 fitting basis with AUXIS XC evaluation option, lower set is GEN-A2 
and BASIS option. See text for details as well as notation for functionals. Eatom and are 
total energies in Hartrees. AE is the total atomization energy in kcal/mol. Re is the equilibrium 
bond length in A. 



Functional 


Eatom 




AE 


Re 


PBE 


-7.460992748 


-14.953949056 


20.06 


2.7236 


PBEMA 


-7.457436406 


-14.946932892 


20.01 


2.7218 


PBEMB 


-7.441633876 


-14.915640294 


20.31 


2.7155 


PBEMC 


-7.452520577 


-14.937217688 


20.19 


2.7196 


PBEMD 


-7.449007723 


-14.930264047 


20.24 


2.7181 


PBEME 


-7.433703689 


-14.899876726 


20.37 


2.7131 


PBE 


-7.460613173 


-14.953310471 


20.13 


2.7304 


PBEMA 


-7.457076546 


-14.946319641 


20.18 


2.7277 


PBEMB 


-7.441303732 


-14.915029100 


20.34 


2.7181 


PBEMC 


-7.452179554 


-14.936623118 


20.25 


2.7244 


PBEMD 


-7.448674591 


-14.929672255 


20.28 


2.7222 


PBEME 


-7.433359029 


-14.899212818 


20.39 


2.7143 


Reference values [20, 21 , 22 , 23] 


-7.47806 


-14.9938 


26 


2.673 



this same size-consistent procedure, namely applying to the separated atoms the modified 
LO constants proper for the value of of the aggregated system (molecule) in question. 
For heteronuclear molecules, especially hydrides, this protocol results in a rather disparate 
enforcement of the LO bound for atoms of substantially different A^, a matter for later 
study and refinement. (We note that the use of the original PBE functional implies the 
most disparate enforcement of all, as it amounts to using the largest N —>■ oo value of A for 
all finite A^.) 

Table IIIII also compares the effect of the two different options in deMon2k for evaluation 
of the XC matrix elements. First is the deMon A2 density fitting basis (deMon2k option 
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AUXIS(2)) and the aforementioned deMon2k option ("AUXIS") for evaluation of XC quan- 
tities using the fitted (auxihary) density on a numerical grid. Second is the richer GEN-A2 
fitting basis and evaluation of the the XC quantities from the density formed straightfor- 
wardly from the KS orbitals also on the numerical grid ("BASIS" option). In principle, the 
latter procedure is the more accurate and is the one we adopted. Nevertheless, the trends 
in the PBEMx series are essentially the same in the less-accurate procedure. 

Our other exploratory test was the O2 molecule, a triplet ground state system. The 
TZVP PEE atomization energy (see the first two lines of Table IIV|) is about three percent 
off from the published result of Ernzerhof and Scuseria jiol], who used the substantially richer 
basis 6311+G(3df,2p). An ACES-II 2J| calculation using another rich basis (aug-cc-PVTZ) 
matched the deMon2k results with that same KS basis and the richer density-fitting basis 
("GEN-A2*" option). These results, in the third and fourth lines of Table HVt calibrate the 
effects of basis set differences. Results for the PBEMx series in the aug-cc-PVTZ/GEN-A2* 
basis sets also are in Table |Vl The relative shifts among the six PEE variants are the same 
irrespective of basis sets, but the richer basis sets make the atomization energies larger and 
bond lengths slightly shorter. Ernzerhof and Scuseria 191] quote the experimental data as 
AE = 118 kcal/mol and = 1.208 A. Thus, all six PBE variants (original plus five new 
give too deep a binding energy at slightly elongated bond lengths. 



TABLE IV: O and O2 (triplet) total energies (Eh), molecular atomization energy (kcal/mol), 
and equilibrium bond length (A) for the ordinary PBE XC functional as calculated in deMon2k, 
TZVP basis ("deMon-TZVP"); GaussianOS, 6311+G(3df,2p) basis (Ref. Q "ES"); ACES-II, aug- 



cc-PVTZ (Ref. 



25| "AP"); deMon2k, aug-cc-PVTZ basis ( "deMon-aug" ) . 



Calc. 



E-, 



tot,0 



E, 



tot,02 



AE 



deMon-TZVP 

ES 

AP 

deMon-aug 



-75.00612438 

-75.00773627 
- 75.00781596 



-150.23282532 

-150.24372619 
-150.24387554 



138.4 
143 
143.2 
143.2 



1.23491 
1.217 
1.21996 
1.22008 



The preceding discussion makes clear that systematic comparison of the five PBE variants 
generally does not require a fully saturated basis set. Rare gas dimers, discussed below, are 
an exception. Therefore, except for rare gas dimers, we adopted the following protocol: 
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TABLE V: As in Table |TTI] but for triplet O2 and for two different basis sets TZVP/GEN-A2 
(upper set), aug-cc-PVTZ/GEN-A2* (lower set). 



Functional 


-^atom 


E02 


AE 


Re 


PBE 


-75.00612438 


-150.23282532 


138.4 


1.2349 


PBEMA 


-74.99689680 


-150.21624772 


139.6 


1.2343 


PBEMB 


-74.95432779 


-150.13915494 


144.6 


1.2315 


PBEMC 


-74.99835094 


-150.21886346 


139.4 


1.2344 


PBEMD 


-74.98900193 


-150.20202494 


140.6 


1.2338 


PBEME 


-74.94639261 


-150.12467907 


145.5 


1.2310 


PBE 


-75.00781596 


-150.24387554 


143.2 


1.2201 


PBEMA 


-74.99853463 


-150.22722985 


144.4 


1.2195 


PBEMB 


-74.95571887 


-150.14986185 


149.6 


1.2167 


PBEMC 


-74.99998015 


-150.22985618 


144.3 


1.2196 


PBEMD 


-74.99059055 


-150.21295118 


145.4 


1.2189 


PBEME 


-74.94774275 


-150.13534295 


150.5 


1.2162 


Reference values [20, 21, 22, 23] 


-75.0674 


-150.2770 


118 


1.208 



(i) use TZVP for the KS basis; (ii) Use GEN-A2 or GEN-A3 algorithms to generate the 
fitting-function basis (and thereby minimize the effect of the variational Coulomb fitting, 
which gives a lower bound to the Coulomb repulsion that can be deceptive with a poorly 
chosen fitting basis); (iii) use the deMon2k BASIS option for evaluation of the XC matrix 
elements. This protocol combines a reasonably rich KS basis with an abundance of caution 
in treating the XC quantities. 

For the rare gas dimers, test calculations on Ar2 with both a DZVP and a 6- 
311++G(3df,3pd) KS basis set demonstrated that these do not reproduce known, large-basis 



PBE results for this dimer 



26 



27l |. Since those two calculations were completely indepen- 



dent and gave essentially identical values, AE = 0.138 kcal/mol. Re = 4.00 A, it is essential 



to reproduce them. Therefore, we shifted to the aug-cc-pVTZ KS basis [28|, as used by Zhao 
and Truhlar, and the deMon2k GEN- A3 fitting function basis. This combination gives the 
same PBE results as the foregoing two references. We treated Ne2 with the corresponding 
aug-cc-pVTZ KS basis and GEN-A3 fitting function basis. 
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IV. RESULTS AND DISCUSSION 



For PBE and PBEMx, Table I VI I gives the atomization energies for eighteen hght 
molecules, while IVIII gives bond lengths and bond angles for those same systems. (As a 
technical aside, note that NiH is a difficult system to treat.) Absolute relative errors in the 
atomization energy are shown in Figure [2] and the corresponding bond length data are in 
Figure [31 

Several features stand out from these results. With a few exceptions, the general pattern is 
that both atomization energies and bond lengths are remarkably insensitive to changes in the 
enforcement of the LO bound. This outcome is consistent with what one might have intuited 
from Figure m The A(A^) which we are justified in using (in the sense that our interpolation 
respects the known constraints) is substantially larger than the A values imputed for actual 
molecules. So one might conclude that a more refined way of implementing the LO bound 
in a GGA is needed. 

However, two systems, Ne2 and Ar2, are notably sensitive to the value of A in the at- 
omization energies. A coherent interpretation by classes of molecules is possible: unlike 
the other atoms, Ne and Ar are closed shell systems. Thus, we may suspect that the well- 
known peculiarities of closed shell interactions are the source of the distinct behavior. For 
the equilibrium bond lengths Ne2 and Ar2 still stand out from all other systems by being 
most sensitive to changes in A. (See Fig. [31) Interestingly, the Be2 bond-length variation 
is at odds with all the other open-shell systems. The behavior of Be2 does not seem to be 
traceable to being from the rapidly varying part of the A(A^) function, since Li2 is in that 
region also and it is insensitive in both bond length and atomization energy. 

To display the effects of imposition of the size-consistent molecular values on the LO 
bounds in atoms, we also calculated some isolated atom total energies at their intrinsic 
N = Z values. For H, C, N, O, and F, Table [VnTl displays the results for the intrinsic value 
versus the results for the highest molecule in which each element was used in the present 
study. As would be expected from the interpolation, Eq. ((51), and the constraints on which 
it is based, the total energy of the H atom exhibits the largest percentage variation between 
intrinsic and molecular values for for the three iV-dependent variants (PBEMC, PBEMD, 
PBEME). 
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FIG. 2: (Top panel): Absolute relative errors in atomization energies for all 18 molecules for 
original PBE and the five variants. (Middle panel): As in the upper panel but with the worst two 
cases (Be2, Ne2) removed to allow a finer scale. (Bottom panel): As in the middle panel but with 
the worst six cases of that panel (F2, Ar2, NiH, Li2, O2, B2) removed to allow a finer scale. 
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FIG. 3: (Upper panel): Absolute relative errors in bond lengths for all 18 molecules for original 
PBE and the five variants. (Lower panel): As in the upper panel but with the two most sensitive 
cases (Ar2, Ne2) removed to allow a finer scale. 

V. CONCLUDING REMARKS 

Our results show that PBE is rather insensitive to changes in A for atoms and covalently 
and ionically bound small molecules. Overall, a reduced, and thus, in principle, better, 
value of A produces slightly worsened energies and slightly improved bond lengths. This 
insensitivity explains why PBE can be successful even though it uses the Aoo value even 
for small A^. In this sense, the present study provides additional insight into the success of 
PBE for small systems. On the other hand, a suitably designed, constraint-based functional 
should give improved results when the constraints it incorporates are sharpened. The failure 
of PBE to meet this expectation must be considered a limitation of the PBE functional form. 

In the case of the closed shell systems, we find a more pronounced A dependence than in 
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TABLE VI: Comparison of effect various Lieb-Oxford bounds in the PBE X functional upon the 
atomization energies (kcal/mol) of various small molecules. See text for notation about functionals. 



Species 


25+1 


PBE 


PBEMA 


PBEMB 


PBEMC 


PBEMD 


PBEME 


Exp [22 


Li2 


1 


20.13 


20.18 


20.34 


20.25 


20.28 


20.39 


26 


Bea 


1 


9.71 


9.95 


10.94 


10.12 


10.35 


11.28 


2.3 


B2 


3 


76.7 


77.3 


80.1 


77.6 


78.2 


80.8 


71.3 


C2 


1 


142.3 


143.4 


147.9 


143.6 


144.6 


148.9 


146 


N2 


1 


235.5 


236.7 


241.7 


236.7 


237.9 


242.6 


227 


O2 


3 


138.4 


139.6 


144.6 


139.4 


140.6 


145.5 


118 


F2 


1 


51.56 


52.46 


56.30 


52.23 


53.12 


56.90 


38 


Ne2 


1 


0.1279 


0.1438 


0.2236 


0.1385 


0.1549 


0.2369 


0.0839 


Ar2 


1 


0.1377 


0.1631 


0.2998 


0.1462 


0.1732 


0.3122 


0.2846 


HF 


1 


142.645 


143.161 


145.401 


143.356 


143.864 


146.037 


142 


LiH 


1 


54.53 


54.54 


54.57 


54.57 


54.57 


54.54 


58 


OH 


2 


110.361 


110.762 


112.458 


110.976 


111.365 


112.981 


107 


NH 


3 


88.95 


89.18 


90.07 


89.34 


89.55 


90.36 


88 


NiH 


2 


74.05 


75.05 


76.63 


74.86 


75.21 


76.79 


58.8 


H2O 


1 


235.9 


236.8 


240.5 


237.1 


237.9 


241.5 


235 


NH3 


1 


305.686 


306.664 


310.761 


307.029 


307.976 


311.882 


297 


CH4 


1 


427.184 


428.245 


432.736 


428.641 


429.673 


433.983 


420 


CeHg 


1 


1423.74 


1430.39 


1458.48 


1425.84 


1432.45 


1460.30 


1362 



the covalent and ionic systems. Because of the dehcate nature of binding in these systems, 
more detailed investigation would be needed to make conclusive statements. In the spirit 
of the preceding paragraph, it would appear to be more productive to focus on developing 
enhancement factors that are more sensitive to sharpening of constraints. 

An appealing thought is that the insensitivity found here may also have to do with the 
way that the LO bound is implemented in DFT in general. The original LO bound is for the 
Coulomb exchange and correlation energy W^c and does not include the correlation kinetic 
energy, = T — T^ > 0, which contributes to E^-c- As a result, E^^ > W^c and the functional 
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TABLE VII: Comparison of effect of various Lieb-Oxford bounds in the PBE X functional upon 
the bond lengths (A) and bond angles (degrees) of various small molecules. The NH3 bond angle 
is Ohnh See text for notation about functionals. CH4 was done with Td symmetry enforced. 



Species 


PBE 


PBEMA 


PBEMB 


PBEMC 


PBEMD 


PBEME 


Exp [23] 


Li2 


2.7304 


2.7277 


2.7181 


2.7244 


2.7222 


2.7143 


2.673 


Be2 


2.4409 


2.4379 


2.4259 


2.4355 


2.4328 


2.4218 


2.45 


B2 


1.6208 


1.6198 


1.6159 


1.6194 


1.6185 


1.6148 


1.590 


C2 


1.2595 


1.2590 


1.2567 


1.2590 


1.2584 


1.2562 


1.243 


N2 


1.1112 


1.1108 


1.1091 


1.1108 


1.1104 


1.1087 


1.098 


O2 


1.2349 


1.2343 


1.2315 


1.2344 


1.2338 


1.2310 


1.208 


F2 


1.4301 


1.4290 


1.4242 


1.4293 


1.4282 


1.4235 


1.412 


Ne2 


3.0808 


3.0418 


2.8709 


3.0546 


3.0162 


2.8550 


3.091 


Ar2 


3.99907 


3.9124 


3.6964 


3.9469 


3.8929 


3.6841 


3.7565 


HF 


0.9385 


0.9385 


0.9385 


09385 


0.9385 


0.9385 


0.917 


LiH 


1.6065 


1.6058 


1.6029 


1.6039 


1.6033 


1.6012 


1.595 


OH 


0.9899 


0.9899 


0.9898 


0.9899 


0.9899 


0.9897 


0.971 


NH 


1.0549 


1.0549 


1.0548 


1.0549 


1.0549 


1.0548 


1.036 


NiH 


1.4580 


1.4594 


1.4553 


1.4599 


1.4590 


1.4549 


1.477 


H2O, R 


0.9750 


0.9749 


0.9747 


0.9749 


0.9749 


0.9746 


0.959 


H20, e 


104.21 


104.24 


104.39 


104.25 


104.29 


104.43 


103.9 


NH3, R 


1.0252 


1.0252 


1.0245 


1.0250 


1.0249 


1.0244 


1.012 


NH3, e 


106.40 


106.44 


106.62 


106.45 


106.49 


106.67 


106.7 


CH4 


1.0990 


1.0989 


1.0984 


1.0988 


1.0987 


1.0982 


1.086 


CgHe, Rcc 


1.3995 


1.3990 


1.3965 


1.3994 


1.3988 


1.3963 


1.397 


CeHg, RcH 


1.0947 


1.0946 


1.0943 


1.0947 


1.0946 


1.0943 


1.084 



X[n] = Exc[n]/E^ [n] which was evaluated in [8|, |9|] is smaller than the functional 

:= . (12) 

If the effect of Tc were large enough, it might explain at least part of the large difference 
between the values of X[n] and X{N) in Figure [TJ What limited numerical evidence we 
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TABLE VIII: Total energies (Hartree a.u.) for five chemically important atoms from various 
Lieb-Oxford bounds in the PBE X functional. Results for the N-dependent functionals are given 
both for the values of N intrinsic to the specific atom and for the highest molecular used: 42 
for H, 42 for C, 14 for N, 16 for O, 18 for F. 



xc 


H 


C 


N 





F 


PBE 


-0 498147969 


-37 794851185 


-54 530389203 


-75 006124382 

1 KJ • \J\J\J J~ ^^KJKJ ^ 


-99 664580137 

iJ kJ '\J\J^.KJKJ\J J^KJ 1 


PBEMA 


-0.497515476 


-37.787528107 


-54.522145702 


-74.996896796 


-99.654626939 


PBEMB 


-0 494726517 


-37 754129975 


-54 484419695 

KJ • ^r^^^m i \j \j \j kj 


-74 954327792 

1 ^ItiJKJ^XKJ^ 1 1 kJ ^ 


-99.608420218 


PBEMC 


-0.481206885 


-37.777286824 


-54.513676785 


-74.989918352 


-99.649171204 


(intTinsir i 

I ±LL\J±. ±LLiJ±\^ J. ^ 1 












PBEMC 


-0.497950959 


-37.792580522 


-54.. 52240237 


-74.99835094 


-99.65716700 


fhiehest A'") 












PBEMD 


-0 481206885 


-37 769885027 


-54 505322010 

KJ ^X. • KJ \J KJ KJ ^ ^ \J -L. \J 


-74 980521867 

1 t %JKJ\JKJ t-i J-KJKJ I 


-99 639010126 


(intrinsic A^) 












PBEMD 


-0.497316032 


-37.785199634 


-54.51406717 


-74.98900193 


-99.64707222 


(highest A^) 












PBEME 


-0.481206885 


-37.736869483 


-54.467831152 


-74.937949507 


-99.592621018 


(intrinsic A^) 












PBEME 


-0.494539297 


-37.751826019 


-54.47640949 


-74.94639261 


-99.60074422 


(highest A^) 












exact [20] 


-0.5 


-37.8450 


-54.5893 


-75.0674 


-99.7341 



have, however, suggests that Aiy[?T.] is only about 10% larger than \\n\, a very modest shift 
compared to the difference in Figure [H 

The more general point, however, that the LO bound is a constraint on exchange and 
correlation together, seems to be sustained by our findings, in that the PBE form enforces 
the bound purely on exchange. One speculation is that the insensitivity found here is in 
part a consequence of that restricted use of the LO bound. 

Finally, we consider aspects of A^-dependence and chemical classification for enforcement 
of the LO bound. The first insight is that, in retrospect, A^-dependent satisfaction of the 
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LO bound actually arose very early in DFT, before the LO proof. In Slater's Xa model, E^c 
is modeled by scaling E^^^. (From a modern perspective, Xa is a one-parameter XC model 



which gains simplicity at the cost of violating correct sea! 



A^-dependent 



33| and exhibits very clear shell structure 



34 



\[n] of Eq. (jl]) is just 3a/2. With typical values of a |3J,|35[, this gives \[n] = 1.0745 for H 



i ng for C.) The a parameter is 



351]. For Xa, the LO functional 



y(A^ = 1) to 1.0387 for Rn (A^ = 86). Comparison with Fig. [T] shows that these values are 
slightly smaller than the highly accurate empirical values found in Refs. j^, 9]. 

The importance of shell-dependent classification was evident in the modern work of Refs. 
Is] and fl. The numerical results of this study also leave a strong suggestion that such 
classification would be helpful. An advantage of the present classification of A with respect 
to A^ alone is that it can be done in an entirely non-empirical way, as it relies only upon 
exact properties of the function A(A^). (This is a clear distinction from all parameterized 
approaches.) Of course, the choice of interpolating function is not unique, but the fact 
that the upper limit on A depends on A^ is completely general. So are all the properties 
of A(A^) used in the construction of our interpolation. What this means is that whatever 
shell-dependent classification might be invented, it must somehow be an addition to (or 
incorporate) classification by particle number, not supplant it. Because that classification 
will have to avoid size-inconsistency, we suspect that the formulation will require additional 
insight, including additional constraints. 
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